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La presente invention concerne une methode pour determiner 1'indice de resistivite 
RI en fonction de la saturation en eau, de certaines roches de porosite complexe. 

La methode selon Tinvention trouve des applications notamment dans le domaine 
5 de l'exploitation de gisements d'hydrocarbures. 

devaluation des reservoirs carbonates est une tache particulierement difficile pour 
les petrophysiciens, et Ton manque toujours de connaissances precises concernant les 
proprietes de transport au sein de ces milieux poreux. Par rapport aux roches siliclastiques, 
les carbonates peuvent etre plus simples sur le plan mineralogique, mais ils sont 

10 incomparablement plus complexes en termes de structure poreuse et de proprietes de 
surface. L'origine largement biologique des sediments, alliee a differents processus 
diagenetiques, donne lieu a des structures poreuses complexes qui peuvent etre 
extremement differentes d'un reservoir a l'autre. Pour de nombreux systemes carbonates, 
les etalonnages-de donnees de resistivite realises en laboratoire sont en contradiction avec 

L5 les observations faites sur le terrain (production anhydre, pression capillaire) et avec les 
mesures directes de saturation en eau effectuees sur des carottes preservees. 

Etat de la technique 

L'etat de la technique est d6fini notamment par les publication suivantes : 

Bouvier L. et al, Reconciliation of Log and Laboratory Derived Irreductible Water 
20 Saturation in a Double Porosity reservoir, Advances in Core Evaluation, edited by 

Worthington and Longeron, Gordon and Breach Science Publishers ; 

- Dixon J.R. et al (1990), The Effect of Bimodal pore Size Distribution on Electrical 
Properties of some Middle Eastern Limestone, Soc. Petr. Eng. 20601, 7 th SPE Middle 
Eastern Oil Show, Bahrain, pp. 743-750 ; 

25 - Fleury M. (1998), 'TRIM : a Fast Resistivity Index Measurement Method", 
Proceedings of the International Symposium of the SZociety of the Core Analysts, Den 
Hague ; 
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Fleury M. et al (2000) "Frequecncy Effect on Resistivity Index Curves Using a New 
Method", Proceedings of the 41 st Annual SPWLA Symposium, Dallas ; 

Moore C.H., (2001) "Carbonate reservoirs, Porosity Evaluation and Diagenesis in a 
Sequence Strati graphic Framework", Developments in Sedimentology 55, Elsevier 
5 Editions ; 

- Petricola M.J.C. et al, (1995), "Effect of Microporosity in Carbonates: Introduction of a 
Versatile Saturation Equation", Soc. Petr. Eng. 29841, SPE Middle Eastern Oil Show, 
Bahrain, pp.607-615 ; 

- Sen P.N. et al (1997), "Resistivity of Partially Saturated Rocks with Microporosity", 
10 Geophysics, Vol.62, N°2, pp. 415-425. 

La comprehension et la prediction de l'effet de la structure du reseau poreux de la 
mouillabilite et des proprietes electriques de roches carbonatees, constituent un veritable 
defi scientifique, tant sur le plan theorique qu' experimental. La correcte evaluation de ces 
parametres a en effet une incidence majeure sur 1'estimation du petrole en place, 
15 notamment pour les gisements petroliers geants du Moyen-Orient, car 1'ecart par rapport 
aux valeurs standard des exposants d'Archie m et n par rapport a la valeur 2 est si 
important que 1'estimation de la saturation en eau peut varier de plus de 20 %. 

Differentes observations experimentales ont montre l'existence de populations de 
pores distinctes: micropores, macropores et mesopores, presentant des degres de 

20 coexistence divers. En general, la courbe d'indice de resistivite RI(Sw)=Rt(Sw)/Ro, ou Rt 
est la resistivite de la roche a une saturation en eau Sw, et Ro la resistivite pour Sw=l, ne 
peut pas etre decrite par une loi de puissance (seconde loi d'Archie RI=Sw" D ) et, cad n est 
fonction de la saturation elle-meme. La microporosite peut agir comme un chemin 
parallele pour le courant, ce qui entraine une diminution des valeurs de n et, de ce fait, une 

25 insensibilite graduelle de la resistivite a la saturation, cornme cela a ete observe sur des 
sables argileux. La microporosite peut etre aussi la cause des faibles valeurs de n mesurees 
(typiquement 1,45). On a aussi observe que n peut augmenter considerate ement dans 
certaines conditions, et qu'il existerait un lien entre la courbe RI(Sw) et la courbe de 
pression capillaire. L' augmentation des valeurs de n est egalement un effet connu de la 

30 mouillabilite, qui tend a favoriser la discontinuity de la phase aqueuse et ainsi a augmenter 
la resistivite par rapport a la mouillabilite a Teau. L'effet de la mouillabilite peut entramer 
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soit une augmentation soudaine de n, soit une valeur n elevee sans discontinuity. Une 
confusion est de ce fait possible. 

Les courbes experimentales RI(Sw) qui ont ete etablies a la suite de divers travaux 
anterieurs deja cites et de nos propres observations, peuvent presenter (Fig.l) quatre 
formes distinctes qui ne satisfont pas toujours aux lois d' Archie : 

- type I : peut etre typique de carbonates proven ant de la formation Thamama, 

- . type II : se redressant a une saturation intermediaire et s'aplatissant a faible saturation 
(la presente etude), 

type m : pente unique a faible saturation, extrapolation a Sw=l au-dessus de In=l, et 

type IV : typique de systemes mouillables a Fhuile, valeurs elevees de n pouvant 
augmenter encore a faible saturation. Ceci est egalement valable pour les clastes. 

On voit done que lorsqu'on se trouve en presence de roches a la porosite complexe 
qui ne satisfont pas aux lois d' Archie, un grand nombre de mesures couteuses est 
necessaire pour tenir compte de la variability de la structure poreuse. 

La methode selon I'invention 

La methode selon l'invention a pour objet de determiner les variations de Tindice de 
resistivite (RI) d'une farnille d f echantillons de roche a structure poreuse complexe en 
fonction de la saturation en eau (Sw), en presence d'un fluide non conducteur. Elle 
comporte les etapes suivantes : 

pour chaque echantillon de la farnille comprenant au moins un premier et un deuxieme 
reseaux de pores, on determine la fraction volumique occupee par chaque reseau de 
pores, par application aux differents echantillons d'une technique de relaxometrie de 
type RMN ; 

pour chaque echantillon de la farnille, on mesure par injection de mercure la 
distribution des seuils de pores dans les differents reseaux de pores ; 
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on determine experimentalement sur un echantillon au moins de la famille servant de 
reference, les valeurs de coefficients reliant la variation de sa resistivite electrique en 
fonction de sa saturation en eau ; et 

- on determine l'indice de resistivite (RI) de tous les echantillons de la famille en se 
basant sur la variation de parametres decrivant 1'agencement du reseau poreux et en 
utilisant les valeurs des coefficients mesures sur l'echantillon de reference. 

Suivant un premier mode de mise en oeuvre convenant pour les echantillons 
comportant deux reseaux de pores, on determine a partir de l'echantillon de reference, les 
valeurs des coefficients reliant la conductivity totale de l'echantillon aux conductivity des 
premier et deuxieme reseaux de pores et aux saturations respectives en eau des deux 
reseaux, l'indice de resistivite etant calcule a partir des fractions volumiques respectives 
des deux reseaux de pores et de la valeur de la saturation moyenne a partir de laquelle le 
reseau ayant les plus petits pores, est envahi par le fluide non conducteur. 

Suivant un deuxieme mode de mise en oeuvre convenant pour les echantillons 
comportant un troisieme reseau de pores, on determine a partir de Techantillon de 
reference, les valeurs des coefficients reliant la conductivity totale de l'echantillon aux 
conductivity des premier et deuxieme reseaux de pores et aux saturations respectives en 
eau des deux premiers reseaux, l'indice de resistivite etant calcule a partir des fractions 
volumiques respectives des trois reseaux de pores, de la valeur de la saturation moyenne a 
partir de laquelle le reseau parmi les deux premiers reseaux de pores qui a les plus petits 
pores est envahi par le fluide non conducteur, et de la valeur de la saturation moyenne a 
partir de laquelle le reseau ayant les plus grands pores parmi les deux premiers reseaux de 
pores, est envahi par le fluide non conducteur. 

Les variations de l'indice de resistivite obtenues par application de la methode, 
s ? ajustent au mieux aux courbes experiment ales dont les exemples sont illustres a la Fig.l. 
La methode s'appuie sur une technique experimentale puissante permettant de mesurer une 
courbe d'indice de resistivite continue au lieu de la technique classique qui decrit la courbe 
RI seulement par un nombre limite de points a l'equilibre capillaire. 
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Presentation succincte des figures 

Les caracteristiques et avantages de la methode et du dispositif selon l'invention, 
apparaTtront plus clairement a la lecture de la description ci-apres en se referant aux dessins 
annexes oii : 

la figure 1 montre une representation schematique en coordonnees logarithmiques de 
differents types de courbes RI en fonction de la saturation en eau Sw observees 
experimentalement; les courbes I a III sont dus a la structure poreuse, la courbe IV est 
typique des effets de la mouillabilite et n'est pas specifique aux carbonates ; 

la figure 2 montre une distribution des tailles de pores de trois carbonates etudies : 
(carbonate de gisement, calcaire d'Estaillade et calcaire de Brauvillier) obtenue par 
RMN ; les chiffres indiquent la fraction volumique de pore a gauche du segment 
vertical ; 

- les figures 3a, 3b montrent la sensibilite du modele a deux reseaux de pores (DPC) a la 
saturation Sc a laquelle les micropores sont envahis par l'huile (3a) et a a, rapport de 
conductivity initiale des deux populations de pores (3b) ; 

les figures 4a, 4b montrent 1'application du modele DPC a l'echantillon RC a la 
temperature ambiante (4a) et dans les conditions de gisement avec de I'huile morte 
(4b); la courbe experimentale est en traits fins et le modele en trait gras; les parametres 
ni, n 2 et a sont ajustes, fi et Sc sont mesures ; 

la figure 5 montre un scenario d 'invasion dans le modele a trois reseaux de pores 
(TPC) ; 

- la figure 6 montre un agencement electrique modelisant les trois reseaux dans le 
modele TPC ; 

- la figure 7 montre une courbe RI mesuree (trait fin) et modele (trait gras) pour 
l'echantillon EL (calcaire d'Estaillade) ; 

la figure 8 montre une courbe RI mesuree (trait fin) et modele (trait gras) pour 
l'echantillon BL (calcaire de Brauvillier) ; et 
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la figure 9 montre une sequence de calibrage simplifie pour les diagraphies de 
resistivite ; pour les roches de type non-Archie, Sw(RI) est necessaire ; les donnees de 
laboratoire fournissent generalement la relation RI(Sw). 

Description detaillee 

On dispose d'un certain nombre d'echantillons d'une meme famille et Ton veut 
calculer Tindice de resistivite de chacun d'eux. La methode que Ton va decrire permet 
d'obtenir ces indices, essentiellement a partir : 

de mesures experimentales effectuees sur chaque echantillon de la famille, donnant le 
nombre de reseaux de pores de tailles differentes et la fraction volumique occupee par 
chaque reseau de pores ainsi que la distribution des seuils de pores pour chacun d'eux ; 
et 

d'une mesure experimentale faite sur un seul echantillon de la famille, servant de 
reference, de la variation de sa resistivite electrique en fonction de sa saturation en eau 
(Sw). 

Nomenclature 

Cto = conductivity totale a Sw=l 

Q = conductivity de la i-ieme population a Swi=l 

fi = fraction volumique de pores RMN de la i-ieme population, i— 1,2,3 

RI = indice de resistivite, Rt/RoSw = saturation moyenne en eau 

Swi = saturation en eau de la i-ieme population de pores 

Sc = saturation moyenne a laquelle les micropores sont envahis 

Sm = saturation moyenne a laquelle les macropores sont envahis 

a h2 = rapport de conductivity Ci/C 2 et C\/C 3 respectivement. 
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Donnees experimentales sur echantillons 

Mesures de temps de relaxation RMN 

Dans un premier temps, on determine par une methode de relaxometrie de type 
RMN bien connue des specialistes, les differents reseaux de pores a I'interieur des 
echantillons de roche et la fraction volumetrique occupee par chaque reseau. Un exemple 
d'une telle methode de relaxometrie d RMN est decrit par exemple dans la demande de 
brevet fran?ais EN 02/... du demandeur. 

Nous presentons les mesures realisees sur trois echantillons de carbonate 
selectionnes. En leur appliquant une methode de relaxometrie de type RMN bien connue 
des specialistes, on peut voir que tous les echantillons se caracterisent par une distribution 
des tailles de pores double ou triple. Le premier est un echantillon de carbonate de 
gisement provenant du Moyen-Orient classe parmi les "packstones". II a ete nettoye au 
moyen de differents solvants a temperature elevee avant de proceder aux mesures. A 
conditions ambiantes, l'experience a ete effectuee avec de l'huile raffinee et de la saumure 
de reservoir synthetique. L' experience realisee dans des conditions de gisement a ete 
conduite avec de l'huile de gisement degazee morte filtree a la pression de reservoir. Apres 
nettoyage, la mouillabilite de Techantillon se caracterise par une mouillabilite moderee a 
1'eau et, apres vieillissement, par une forte mouillabilite a l'huile. Les deux autres 
echantillons sont des carbonates d'affleurement (mouillables a 1'eau). La porosite du 
calcaire de Brauvillier (BL) est essentiellement intergranulaire et due au cortex de 1'oolite. 
La porosite du calcaire d'Estaillade est a la fois inter et intragranulaire. Les populations de 
pores sont separees par un facteur 10 au moins. 

Mesures d'indice de resistivite 

On mesure ensuite Tindice de resistivite des echantillons. On peut utiliser 
avantageusement la methode de mesure rapide de Findice de resistivite (dite FRIM) telle 
que decrite dans les brevets EP 974 839 du demandeur. 

Suivant cette methode, on realise un deplacement force huile-eau sur une petite 
carotte de 2,5 cm de long et 4 cm de diametre par exemple. Ce deplacement est tres proche 
d'un procede de deplacement a paroi semipermeable, excepte que l'equilibre capillaire 
n'est pas necessaire. A conditions ambiantes, seuls deux ou trois paliers de pression sont 
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utilises. Aux conditions de gisement, lors de l'utilisation d'huile morte ou brute, on fait 
durer chaque palier de pression le temps necessaire pour obtenir la stabilisation de la 
mouillabilite a une saturation donnee. Le drainage s'effectue ainsi en approximativement 
deux a trois semaines, ce qui correspond de maniere typique a la cinetique des processus 
5 chimiques lies au vieillissement. La simplicity de cette methode tient au fait qu'il suffit 
d'enregistrer et de reporter en temps reel la saturation et la resistance moyennes pour 
obtenir une courbe d'indice de resistivite continue exempte de tout artefact. Le point-cle 
reside dans le fait que la geometrie radiale des electrodes permet d' examiner la totalite du 
volume de Techantillon et de compenser le profil de saturation non-uniforme qui apparait 
10 en l'absence d'equilibre capillaire. 

A conditions ambiantes, nous avons utilises la cellule a electrodes radiales decrite 
par exemple dans le brevet EP-A-974 839 deja cite. Les mesures d'impedance complexe 
ont ete effectuees a une frequence fixe de 1 kHz et la partie reelle a ete extraite pour 
calculer l'indice de resistivite. La pression capillaire imposee la plus elevee etait de 12 bars 
15 (pour une tension interfaciale y=35 mN/m ; aux conditions de gisement, avec de Thuile 
brute, la pression capillaire maximale a ete reduite proportionnellement a la tension 
interfaciale). 

Modeles de conductivity 

Modele de conductivite de porosite double (DPC) 

20 L'objectif est d'expliquer le flechissement des courbes RI a faible saturation (type 

I, figure 1). Le modele de double porosite considere ici est pour l'essentiel tres proche de 
ceux proposes pour les sables argileux, ou les argiles presentes a la surface des pores 
constituent un chemin parallele pour le courant. Pour les carbonates, nous admettons 
Texistence de deux reseaux de pores presentant des conductivites electriques paralleles. 

25 Les deux principaux ingredients dans notre modele sont la description de l'invasion des 
reseaux de pores en cours de drainage, et la description de l'agencement electrique des 
deux differentes populations de pores. 

Le premier reseau 1 (macropores par exemple) represente la majeure partie du 
volume de pores, tandis que le deuxieme reseau 2 (micropores) n'en represente qu'une 
30 petite fraction, qui n'est pas necessairement superieure aux seuils de percolation. Dans un 



1er depot 

9 



premier temps, nous considerons la saturation de chaque reseau, Swi et Sw 2 , qui sont 
reliees a la saturation moyenne mesuree Sw selon : 

Sw = fjSwj + f 2 Sw 2 ou fi -h f 2 = 1 (1). 

fi et f 2 represented la fraction , volumique de pore de chaque population. Ces 
fractions ont ete evaluees au moyen de la relaxometrie RMN. Dans un deuxieme temps, on 
admet que les reseaux sont envahis par de l'huile a differentes pressions capillaires ; l'huile 
accede aux pores de petite taille a une pression superieure a celle observee pour les pores 
de grande taille. A partir de la courbe de pression capillaire etablie precedemment, cette 
pression correspond a une saturation moyenne Sc qui peut etre deduite des courbes 
obtenues par injection de mercure (indiquant la distribution des seuils de pores). Swi peut 
etre exprime en fonction de Sw a saturation Sw elev6e : 

S W] = — + f x ~ 1 , Sw 2 = 1 pour Sw > Sc (2) 
J\ 

En dessous de Sc, les relations Swi=f(Sw) et Sw 2 =f(Sw) necessitent d'autres 
hypotheses. On admettra (i) une relation lineaire et (ii) que Sw^O, Sw 2 ^0 quand 

On en deduit que : 

Sw. =^ Sw, 

Sw 

Sw 2 = pour Sw < Sc (3) 

Sc 

Nous considerons maintenant la conductivite de chaque reseau. Dans le cas d'une 
saturation initiale en saumure, la conductivite totale Cto pour les deux reseaux en parallele 
sera : 

Ct 0 = Ci + C 2 = Ci (1 + a) ou C 2 = a C, (4). 

Le pai-ametre a est le rapport de conductivite des deux reseaux satures a 100 %. A 
partir de la premiere loi d'Archie, on suppose que Ci et C 2 sont relies a la fraction 
volumique de pore de chaque population et que, de ce fait, a est de l'ordre du rapport f 2 /f,. 
Lorsque les deux reseaux sont envahis par l'huile, on admet que chaque conductivite est 
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liee a la saturation par une loi de puissance (comme dans la seconde Joi d'Archie). La 
conductivity totale dans les deux domaines de saturation est la suivante : 

Ct = Sw, nl Ci + C 2 pourSw>Sc (5) 

Ct = Sw! nl d + Sw 2 * 2 C 2 pourSw<Sc (6). 

Au moyen des equations 4, 5 et 6, on trouve l'indice de resistivite RI suivant : 

RI = % = Sw; nl 1 + a : pour Sw > Sc (7) 
Ct l + aSw; nl 

RI = 5 V ! n2 P° ur s ™ * Sc (8) 

Lorsqu'une population de pores domine (le reseau 1 dans le cas present), les 
fonctions RI sont regies par les proprietes de resistivite (m) de cette population. Les 
equations 7 et 8 sont essentiellement similaires aux formules utilisees pour les sables 
argileux, excepte que nous introduisons un second exposant n 2 qui caracterise le second 
reseau. 

Dans le modele DPC, il y a au total 4 parametres ni, n 2 , a et Sc. Pour une courbe 
experimentale donnee, Sc est mesure separement lors d'une experience d'injection de 
mercure-alors que les autres parametres sont ajustes. On observe toutefois un intervalle de 
variation de a autour de f 2 /fi pour lequel il existe une explication physique. A saturation 
Sw elevee, la pente a Fechelle bilogarithmique de RI(Sw) est - n 1} et a faible saturation, la 
pente est - n 2 . Dans une certaine mesure, Sc et a se compensent l'un 1' autre (figure 3), 
mais a est le parametre le plus sensible qui controle la valeur finale de RI. On notera que le 
cas Sc=0,05 presente sur la figure 3 correspond a une situation ou le deuxieme reseau n'est 
pas envahi par Thuile, ce qui donne une asymptote horizontals 

Les parametres du modele ont ete ajustes en vue de correspondre a la courbe 
experimentale RI mesuree sur l'echantillon RC. La valeur de Sc=0,4 a ete deduite du coude 
de la courbe de pression capillaire obtenue par injection de mercure et la fraction 
volumique f,=0,88 de macropores (reseau 1) a ete deduite par RMN (figure 2). Malgre la 
faible permeabilite de Techantillon, la saturation la plus faible obtenue est tres basse (4 %), 
ce qui permet une bonne determination des parametres du modele aux conditions 
ambiantes. La pente initiale de la courbe RI est ni=l,71 ; elle caracterise le reseau 1. Le 



1er depot 

11 



reseau 2 est tres peu sensible a la saturation (n 2 =0,25) et ne se comporte pas comme un 
reseau standard. Toutefois, le rapport de conductivite initiale a=0,054 du reseau 1 a 2 est 
de Tordre de grandeur de f 2 /fi=0,136. Aux conditions de gisement (conditions de 
mouillabilite a Fhuile), on observe une forte augmentation de m, mais la courbe RI n'est 
toujours pas lineaire a Techelle bilogarithmique. Les caracteristiques n 2 =0,78 et a=0,l 1 du 
deuxieme reseau sont aussi legerement modifiees, mais la precision sur ces parametres est 
plus faible qu'avec les conditions ambiantes car la saturation atteinte (a la meme pression 
capillaire) est tres superieure et proche de Sc. 

Formule approchee dans le cas du modele DPC 

Nous proposons une formule approchee reduisant les relations 7 et 8 a une formule 
unique, qui est valide pour les courbes de type I. Nous notons d'abord que : 

Sw 2 n2 Swf nl * Sc" n2 Sw n2 ' Dl « Sw n2 " nl 

avec Swi s Sw et Sc" n2 « 1 (15) 

ce qui donne f 2 «f i . Comme n 2 « ni, les equations 7 et 8 peuvent etre approchees par : 



RI = Sw~ nl -——i — r (1 6). 

l+CSw n2 ~ nl 



Comme le montre la figure 10, Futilisation des 3 parametres ni, n 2 et C suffit pour 
decrire avec precision les donnees au moyen de l'equation 16. Ici, la signification de C est 
approximativement la meme que a et peut dependre de la temperature. La formulation 
donnee dans l'equation 16 peut etre utilisee pour decrire des experiences standard realisees 
a Tequilibre capillaire pour lesquelles on ne dispose que d'un nombre limite de points. 

Modele de conductivite de porosite triple (TPC) 

L'objectif est d'expliquer le redressement des courbes RI a saturation intermediaire 
ou elevee, et leur inflexion a faible saturation a l'echelle bilogarithmique (types II et III, 
figure 1). Nous nous servons ici de Fidee generale selon laquelle les carbonates complexes 
peuvent avoir trois populations de pores appelees, dans un souci de simplification, micro, 
macro et mesopores (3, 2 et 1 respectivement). Comme dans le modele de porosite double, 
on considere la saturation des 3 populations : 
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Sw = f!Sw] + f 2 Sw y +f 3 Sw 3 ou f 1 +f 2 + f 3 =l (9) 

La encore, on admet que l'invasion de ces populations par l'huile au cours du 
drainage est sequentielle. Si le reseau 1 est envahi le premier, on definit une saturation 
moyenne Sm a laquelle le reseau 2 est envahi : 



Sw 2 =1 Sw 3 = 1 pour Sw > Sm (10) 

En dessous de Sm, il est possible d'imaginer de nombreux scenarios. En general, on 
admet des relations lineaires pour les fonctions Swi(Sw) et Sw 2 (Sw). Selon un scenario 
possible, on admettra que Sw a 0, Sw 2 -» 0 lorsque Sw — > Sc. Sc est la saturation a 
laquelle la microporosite est envahie. H s'ensuit que : 



Ce scenario d' invasion est resume sur la figure 5. Typiquement, on pense a une 
situation ou f3«fi et fi-f 2} et la microporosite (reseau 3) est envahie a une pression 
beaucoup trop elevee pour pouvoir etre observee au cours de Texperience. 

Considerons maintenant la conductivite de ces populations. On admet que les 
reseaux 1 et 2 sont en serie, tandis que le reseau 3 est en parallele comme le montre la 
figure 6. La disposition en serie peut sembler en contradiction avec une invasion 
sequentielle par l'huile (qui est un mecanisme essentiellement parallele), mais une telle 
situation est possible. A Sw=l, la conductivite Ct 0 du systeme represente sur la figure 6 est 
la suivante : 




Sw-Sc 
Sm. - Sc ' 



Sw 3 = 1 pour Sc< Sw< Sm 



(ID 



Ct 0 = (Cf 1 + C^ 1 )* 1 + C 3 = C, [(1 + ccf + a 2 ] 



ou C 2 = a] Ci,C 3 = a 2 Ci 



(12). 



Dans les deux domaines de saturation, on obtient : 



Ct = Cj [(Swf nl + af 1 )' 1 + a 2 ] pour Sw > Sm 



(13) 



Ct = d [Swi nl + Sw 2 " q2 af 1 )" 1 + a 2 ] pour Sc < Sw < Sm (14). 
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L'indice de resistivite peut etre calcule a partir des equations 12 a 14. Pour les 
valeurs mesurees de fi, f 2 et f 3 , il convient d'ajuster ni, n 2 , otj et oc 2 aux donnees 
experimentales. Tout comme dans le modele DPC, on s'attend a trouver a } de l'ordre de 
f 2 /fi et a 2 de l'ordre de f^fi. 

Ce modele a ete teste avec des courbes d'indice de resistivite mesurees sur des 
echantillons EL et BL. Sur l'echantillon EL (figure 7), le redressement de la courbe a 
saturation elevee est reproduit qualitativement a saturation elevee et l'aplatissement a 
faible saturation est bien reproduit. On observe une discontinuity a Sm (deduite du coude 
de la courbe de pression capillaire obtenue par injection de mercure) en raison d'une 
modification soudaine de la relation Swj(Sw) (telle que celle montree sur la figure 5). 
Apres ajustage, les deux populations de pores dominantes presentent des caracteristiques 
similaires (n 1 =n 2 =l,5). La troisieme population, qui se comporte comme un circuit 
parallele, n'est pas envahie par Fhuile et elle est seulement caracterisee par le rapport de 
conductivity a 2 . Sur Techantillon BL (figure 8), la courbe RI peut etre expliquee par un 
mecanisme similaire de resistance en serie, mais avec une valeur Sm tres elevee. Nous 
observons enfin que les deux reseaux possedent des valeurs n tres differentes et que le 
reseau 2 est beaucoup moins conducteur que le reseau 1 (a!=7,l). 
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REVENDICATIONS 

1) Methode pour determiner les variations de l'indice de resistivite (RI) d'une 
farnille d'echantillons de roche a structure poreuse complexe en fonction de la saturation en 
eau (Sw), en presence d'un fluide non conducteur, comportant les etapes suivantes : 

5 - pour chaque echantillon de la farnille comprenant au moins un premier et un deuxieme 
reseau de pores, on determine la fraction volumique (f is f 2 ) occupee par chaque reseau 
de pores, par application aux differents echantillons d'une technique de relaxometrie de 
type RMN ; 

pour chaque echantillon de la farnille, on mesure par injection de mercure la 
10 distribution des seuils de pores donnant acces aux differents reseaux de pores ; 

on determine experimental ement sur un echantillon au moins de la farnille servant de 
reference, les valeurs de coefficients (n 1? n 2 ) reliant la variation de sa resistivite 
electrique en fonction de sa saturation en eau (Sw) ; et 

- on determine 1'indice de resistivite (RI) de tous les echantillons de la farnille en se 
15 basant sur la variation de parametres decrivant l'agencement du reseau poreux (fj, f 2 , 

Sc, Sm) et en utilisant les valeurs des dits coefficients mesures sur l'echantillon de 
reference. 

2) Methode selon la revendication 1, caracterisee en ce que Ton determine a partir 
de l'echantillon de reference, les valeurs des coefficients (nj, n 2 ) reliant la conductivite 

20 totale (Ct) de l'echantillon aux conductivite (d, C 2 ) des premier et deuxieme reseaux de 
pores et aux saturations respectives en eau (SW h SW 2 ) des deux reseaux, l'indice de 
resistivite etant calcule a partir des fractions volumiques respectives (fj, f 2 ) des deux 
reseaux de pores et de la valeur de la saturation moyenne (Sc) a partir de laquelle le reseau 
ayant les pores les plus petits, est envahi par le fluide non conducteur. 

25 3) Methode selon la revendication 1, caracterisee en ce que pour un echantillon 

comprenant un troisieme reseau de pores, on determine a partir de l'echantillon de 
reference, les valeurs des coefficients (ni, n 2 ) reliant la conductivite totale (Ct) de 
l'echantillon aux conductivite (C l5 C 2 ) des premier et deuxieme reseaux de pores et aux 
saturations respectives en eau (SWi, SW 2 ) des deux premiers reseaux, Findice de resistivite 

30 etant calcule a partir des fractions volumiques respectives (fi, f 2> f 3 ) des trois reseaux de 
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pores, de la valeur de la saturation moyenne (Sc) a partir de laquelle le reseau qui a les plus 
petits pores parmi les deux premiers reseaux de pores, est envahi par le fluide non 
conducteur, et de la valeur (Sm) a partir de laquelle le reseau ayant les plus grands pores, 
parmi les deux premiers reseaux de pores, est envahi par le fluide non conducteur. 
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